Nanoscale freezing of the 2D spin liquid PraGasSiOu 
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In this letter, we report on the single crystal growth and physical characterization of the distorted 
kagome system Pr3Ga5SiOi4. It is found that at zero magnetic field the system shows no magnetic 
order down to 0.035 K and exhibits a T^ behavior for the specific heat at low temperatures, indicative 
of a gapless 2D spin liquid state. Application of an applied field induces nanoscale islands of ordered 
spins, with a concomitant reduction of the T^ specific heat term. This state could be a possible 
ferro-spin nematic ordering stabilized out of an unusual spin liquid state. 
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One of the desired targets of modern materials sci- 
ence is the realization of model magnetic systems to 
test physical laws. In the field of geometrically frus- 
trated magnetism, the Hamiltonians for many lattices 
have been expressed for various two dimensional and 
three dimensional sublattice, such as the triangular lat- 
tice and the pyrochlores, and experimental analogues 
have been well studied[l|. However, the two dimensional 
kagome lattice has proved to be more difficult to synthe- 
size, as many examples found in the literature (such as 
SrGa8_a;Ga4+a;Oi9, or SCGO) have site disorder issues, 
or are difficult to produce in single crystalline form 2]. 
Jarosites, Such as (D30)Fe3(S04)2(OD)6, have a simi- 
lar problem with a significant amount of Fe vacancies 
for many species of this family 3] . Recently, the parat- 
acamite family Zna;Cu4_2;(OH)6Cli2Q], which shows a 
spin liquid state, was reported. These materials not only 
have "structurally perfect" bond angles for the kagome 
lattice - they are also composed of s = 1/2 Cu^+ spins. 
It has been shown that these materials are excellent test- 
ing grounds for the resonating-valence-bond state pre- 
dicted by Anderson which has relevance not only for 
the study of quantum spin liquids, but also for stud- 
ies of high-temperature superconductors and other low- 
dimensional magnetic systems. However, the issue of site 
disorder of the paratacamite still needs to be addressed. 

Many quantum and classical theories for the kagome 
lattice predict ordered ground states (typically with a \/3 
X Vs structure). The exceptions to this trend include 
SCGO, the jarosite (D30)Fe3(S04)2(OD)6, and the 
newly discovered paratacamites Zn2,Cu4_2;(OH)6Cli2. 
All of these compounds exhibit unconventional spin freez- 
ing at low temperatures. The former compounds exhibit, 
for example, a magnetic component to the specific 
heat (rather than linear component as expected for a 
spin glass) 0. The paratacamite shows a sublinear T 



component that is believed to be the result of a Fermi 
surface of quantum excitations appearing at low tem- 
peratures, although this has yet to be confirmed with 
single crystals]^. All of these samples show dynamics on 
slow time scales that can be resolved with high resolution 
inelastic neutron scattering fzj, muon spin relaxation 
or The origin for these slow spin dynamics is 

thought to lie with persistent two dimensional quantum 
fluctuations which stabilize a spin liquid state, even in 
the case of the larger spins with Cr'^^ (S=3/2) and Fe'^+ 
(S = 5/2). 

In this letter, we report the single crystal growth 
and characterization of the distorted kagome system 
Pr3Ga5SiOi4. This material has a two dimensional lat- 
tice of Pr'^+ spins (J = 4) with weak antiferromagnetic 
couphng (6'cw = —2.3 K). Neutron scattering experi- 
ments fail to show any ordering in zero applied magnetic 
flelds down to 35 niK, despite the presence of a broad 
peak in the heat capacity at T = 6.7 K. The presence of 
a component in the heat capacity is correlated with 
low energy two dimensional spin excitations. Applica- 
tion of an external magnetic field perpendicular to the 
two dimensional layers induces a gap in the spin excita- 
tion spectrum, but magnetic Bragg peaks indicative of 
ordering do not appear. There is only diffuse scattering 
present up to H = 9 T which is modeled with the forma- 
tion of nanoscale islands of ordered spins. This sort of 
ordering has been seen in other two dimensional systems 
such as NiGa2S4[l3|, but has only now been observed in 
a kagome lattice with applied fields. This state could be 
a possible ferro-spin nematic ordering stabilized out of 
an unusual spin liquid state. 

A single crystal of Pr3Ga5SiOi4 was grown by the 
traveling-solvent fioating-zone technique. The room- 
temperature structure was determined with an X-ray 
diffractometer equipped with Cu Kqi radiation. X-ray 
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FIG. 1: (a) Distorted kagome lattice of Pr'^+ ions in the 
ab plane of Pr3Ga5SiOi4. (b) Temperature dependences 
of DC susceptibility and its inverse, (c) XRD pattern for 
Pr3Ga5SiOi4 (plus marks) at room temperature. The solid 
curve is the best fit from the Rietveld refinement using Full- 
Prof. The vertical marks indicate the position of Bragg peaks, 
and the bottom curve shows the difference between the ob- 
served and calculated intensities. Inset: Neutron diffraction 
patterns at T = 0.035 K, H = T; T = 0.035 K, H = 9 T; 
and r = 40 K, H = T show no magnetic Bragg peaks. 

Laue diffraction was used to orient the crystal. The 
magnetic-susceptibility measurements were made with a 
DC superconducting interference device (SQUID) mag- 
netometer; the measurements were made on heating after 
cooling in zero field and with applied magnetic field (H = 
1 T) parallel to the longest dimension of the sample. The 
specific heat and AC susceptibility measurements were 
made with a physical property measurement system on 
single crystals with applied fields parallel to c axis. A di- 
lution fridge was used to collect specific heat data below 
0.4 K. Neutron scattering measurements were completed 
at the NIST CHRNS using the Disk Chopper Spectrom- 
eter with a wavelength of 5.0 A. The crystal (total mass 
of 5 g) was aligned in the ab plane with a vertical mag- 
netic field applied in the c direction. A dilution fridge 
was used which had a base temperature of 0.035 K 

Pr3Ga5SiOi4 crystallizes in the trigonal space group 
P321 with lattice parameters a — 8.0661(2) A and c 
= 5.0620(2) A (Fig. 1(c)). The Ft^+ magnetic ions in 
Pr3Ga5SiOi4 are organized in corner sharing triangles in 
well-separated planes perpendicular to the c axis. Within 
each plane, the Pr^+ ions form a distorted kagome lat- 
tice (Fig. 1(a)), which is topologically equivalent to the 
ideal kagome when only the shortest atom bridging in- 
teractions are considered. The nearest Pr-Pr separation 
in the ab plane is 4.2 A (Fig. 1(a)). 

Figure 1(b) shows the temperature dependences of 



FIG. 2: (a) Temperature dependences of Cp/T for 
Pr3Ga5SiOi4 and La3Ga5SiOi4. (b) Temperature depen- 
dences of Cmag/r and the integrated entropy for Pr3Ga5SiOi4 
and the crystal field calculation. Open circles are experimen- 
tal results, solid lines are calculated results, (c) Tempera- 
ture dependences of the specific heat at low temperatures for 
Pr3Ga5SiOi4 with H = T, 6 T, and 9 T. (d) Magnetic contri- 
bution of specific heat after subtracting the Schottky anomaly 
and lattice contribution. 

the DC magnetic susceptibility (xiT)) and its inverse 
(x^^(r)). The susceptibility above 50 K follows the 
Curie- Weiss law: x'^T) = (T - 6'cw)/C. The effec- 
tive moment ficS = 3-2(1) /^^b calculated from the Curie 
constant C is smaller than that of single Pr^+ ion (/ics 
= 3.6 /is), which is common for magnetically frustrated 
systems[l[. The Curie- Weiss temperature, 6'cw = —2.3 
K, indicates that the nearest neighbor interactions are 
antiferromagnetic. No magnetic anomalies were observed 
down to 1.8 K from the DC susceptibility data. The in- 
set of Fig. 1 (c) shows the elastic neutron diffraction 
patterns at different temperatures and magnetic fields. 
No magnetic Bragg peaks or significant magnetic diffuse 
scattering features are found down to 0.035 K. Note that 
this gives a frustration index of / ~| 6'cw I /Tc > 66. 

The specific heat data, Cp (T) , also shows no evidence 
of a phase transition down to 0.1 K in zero field (Fig. 2(a) 
and (c)). The magnetic specific heat Cinag{T)/T af- 
ter substraction of the lattice contribution (Cp(T) of 
La3Ga5SiOi4 with no magnetic ions) exhibits a broad 
peak at Tpoak — 6.7 K. The integration of Cmag/T gives 
an entropy S'mag = 8.5 J/mol-K ( circles in Fig. 2(b)) near 
200 K, which is approaching i?ln3 {R is the gas constant). 
This value suggests that the crystal field scheme of Pr^"*" 
below 200 K could be (i) three singlet levels, (ii) a dou- 
blet ground state and a singlet excitation level, or (iii) a 
singlet ground state and a doublet excitation level. These 
three possibilities have been tested with the experimen- 
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tal data - the best fit is obtained from the three singlet 
level scheme with Ai = 25 K and A2 = 117 K. This is in 
good agreement with earlier studies on the related gar- 
net structure Pr3Ga50i2, which has three singlets with 
an energy separation of Ai = 26 K and A2 = 68 K 11 1 . 
It would be unusual if Pr3Ga5SiOi4, of lower site sym- 
metry, would have doublets if the cubic garnet structure 
is already composed of low temperature singlets. The 
next singlet of Pr3Ga50i2 is at T = 780 K, and this is 
not seen in our specific heat data. The calculated data 
(the line in Fig. 2(b)) fits experimental data well at high 
temperatures; but at low temperatures the experimental 
data shows a more broad peak and larger value of spe- 
cific heat below Tpeak- The origin for this feature is the 
broadening of the crystal field levels due to correlation 
effects, as seen in other spin liquid candidates such as 

Tb2Ti2 07fl2t 

At low temperatures (Fig. 2(c)), the nuclear Schot- 
tky anomaly introduces an upturn in the specific heat, 
which can be fit with a term. After the substrac- 
tion of this anomaly and the lattice contribution, the 
magnetic contribution of the specific heat at zero field 
exhibits power-law behavior at low temperatures (Fig. 
2(d)). The data between 0.1 K and 4 K are well fitted 
by a power law Cmag — AT"', where A is a constant 
and a = 1.98(2). Two dimensional (2D) spin excitations 
would give C ^ T^. This quadratic temperature depen- 
dence without long-range magnetic order indicates the 
presence of gapless linear modes in 2D, similar to other 
2D kagome systems such as SCGO.Q 

In general, a peak in C(T) /T results from a peak in the 
density of states, g{w), defined as U = J dwg{w)n{w)w, 
where U is the internal energy, n{w) is the Bose pop- 
ulation factor, and the integral is taken over the excita- 
tion bandwidth. Usually, the temperature of the C{T)/T 
peak is roughly half the mode energy. Applying this rule 
here, a peak in g{w) is expected at hwo ~ 13 K. In order 
to confirm this prediction, inelastic neutron scattering 
experiments were completed. In Fig. 3(c), the results 
of the integrated intensity scans are shown as a function 
of applied field. Note that at zero field there is a broad 
peak at E = 1.2 meV ^ 13 K, which is consistent with 
the temperature of the specific heat peak. 

In summary, the zero field data on single crystalline 
Pr3Ga5SiOi4 shows the following low temperature prop- 
erties: (i) the absence of long-range magnetic order; (ii) 
the absence of magnetic diffuse scattering; (iii) a de- 
pendence of the specific heat; and (iv) the existence of 
spin excitations related to a highly degenerate state as 
T approaches zero. These properties place strong con- 
straints on possible ground states. Observations (i), (iii), 
and (iv) are consistent with a spin liquid with no con- 
ventional long-range magnetic order, such as observed 
in the triangular lattice NiGa2S4[l3l and hyper- kagome 
Na4lr308'l3'|. However, the absence of magnetic diffuse 
scattering is unusual. One possibility is that the spins 
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FIG. 3: Temperature dependences of the AC susceptibility: 
(a) real part and (b) imaginary part at frequency 1 Hz and 
1000 Hz. (c) Inelastic neutron scattering integrated over Q at 
0.035 K as a function of magnetic field. 



are truly dynamic on the neutron time scale. The AC 
susceptibility (Fig. 3(a, b)) shows no anomaly, nor any 
frequency dependence down to 1.8 K, which also indi- 
cates that the characteristic fluctuation rate is beyond 
the kHz region. The featureless AC susceptibility data 
also excludes the possibility that Pr3Ga5SiOi4 is a spin 
glass at low temperatures. 

The specific heat measured at magnetic fields 
(Fig. 2(d)) shows that the absolute value becomes smaller 
and behavior disappears. These field-dependent be- 
haviors can be related to the inelastic spin excitation 
(Fig. 3(c)), which is clearly suppressed by applied fields. 
The suppression of the density of states of the 2D spin 
excitations in higher fields is clear from the neutron scat- 
tering data. The elastic neutron pattern at 0.035 K with 
H = 9 T (Inset of Fig. 1(c)) shows no extra peak and in- 
tensity compared to the zero field pattern, which means 
no long-range magnetic order is stabilized. However, 
significant diffuse scattering appears near Q = 0.78 A 
(Fig. 4(a)) with H = 9 T. The form of this diffuse scat- 
tering of a broad asymmetric peak (a sharp rise at low 
Q and a slow fall off toward high Q) is characteristic of 
two-dimensional short-range order. The scattering from 
Pr3Ga5SiOi4 can be described analytically by a modified 
Warren function for 2D magnetic correlations 1J| . The 
structure factor P20 , around the peak is expressed by: 
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FIG. 4: (a) DifTuse scattering of PraGagSiOw at 0.035 K 
with different magnetic fields, (b) Difference between patterns 
taken at H = 6 T , 9 T and that taken at H = T. The solid 
lines are fits to Eq. 1. Inset: Schematic plot of magnetic 
clusters with size L. 
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clusters induces a spin gap in the excitation spectrum, 
and subsequently reduces the component of the spe- 
cific heat. 

For the studied geometrically frustrated systems, such 
as SCG0[2], NiGa2S4[iO|], and Na4lr308[l3], one com- 
mon feature is that the behavior of specific heat is in- 
dependent of the applied magnetic field, suggesting that 
the ground state consists of moment free spin clusters. 
In contrast, Pr3Ga5SiOi4 studied here has a behavior 
which is sensitive to magnetic fields with a gap opening in 
the spin excitation spectrum. One of the possible states 
which could correspond to such behavior is the recently 
predicted spin nematic that shoul d ap pear in two dimen- 
sional triangule-based lattices l^, 1^. However, most of 
these models in the literature lack terms in the Hamilto- 
nian to express the crystal field effects of the Pr3+ sites 
(they are modeled on Heisenberg like transition metal 
systems with quenched orbital angular momentum). De- 
spite this, we can rule out the presence of a non-coUinear 
spin nematic since much of the spin excitation spec- 
trum is gapped out within the ah plane. A ferro-nematic 
state, as mentioned by Bhattacharjee et a^flij. in which 
the moment is spatially uniform and the excitations are 
polarized, is a more likely possibility. Future experiments 
are clearly needed to verify or refute this scenario. 
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a = (2V^i/A)(sin0 - sin6'o) (2) 

and with 

/■lO 

F{a)^ exp[-(a;2--a)2]da; (3) 
Jq 

L is the spin-spin correlation length, A is the wave- 
length, if is a scaling constant, m is the multiplicity of 
the reflection, Fhk is the two-dimensional structure factor 
for the spin array, and is the center of the peak. The 
fits of the diffuse scattering to equation (1) (Fig. 4(b)) 
give a correlation length L = 25(1) A and 29(1) A for 
H = 6 T and 9 T, respectively. This correlation length 
corresponds to 6 ~ 7 in-plane kagome lattice spacings. 
The increasing L with increasing field shows that the 
short-range order or the nanoscale magnetic cluster size 
expands with increasing fields. The formation of these 
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